Quantum diffusion with disorder, noise and interaction 
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Disorder, noise and interaction play a crucial role in the transport properties of real systems, but they are typ- 
ically hard to control and study both theoretically and experimentally, especially in the quantum case. Here we 
explore a paradigmatic problem, the diffusion of a wavepacket, by employing ultra-cold atoms in a disordered 
lattice with controlled noise and tunable interaction. The presence of disorder leads to Anderson localization, 
while both interaction and noise tend to suppress localization and restore transport, although with completely 
different mechanisms. When only noise or interaction are present we observe a diffusion dynamics that can 
be explained by existing microscopic models. When noise and interaction are combined, we observe instead a 
complex anomalous diffusion. By combining experimental measurements with numerical simulations, we show 
that such anomalous behavior can be modeled with a generalized diffusion equation, in which the noise- and 
interaction-induced diffusions enter in an additive manner. Our study reveals also a more complex interplay 
between the two diffusion mechanisms in regimes of strong interaction or narrowband noise. 



Disorder, noise and interaction are known to play a funda- 
mental role in the dynamics of quantum systems, but a de- 
tailed understanding of their combined action is still missing. 
The interest in this general problem ranges from electronic 
systems [1], spin glasses [2J, and nanoscale quantum Brow- 
nian motors [3|, to quantum communication (4J [5) and the 
physics of biological complexes llrJ UTUl . Despite its impor- 
tance, there is so far only a very limited theoretical under- 
standing of the interplay of disorder, noise and interaction. 
Moreover, these ingredients are hard to control in experiments 
with natural or artificial systems; one exception are quantum 
optical schemes iflTj - TPfl . where nonlinearities are however 
weak and do not allow to investigate in depth many-body ef- 
fects. 

The prototypical dynamical problem for disordered quan- 
tum systems is the evolution of an initially localized 
wavepacket [15]. To our knowledge this problem has never 
been studied, neither theoretically nor experimentally, under 
the combined effect of noise and interaction. It is well known 
that in one spatial dimension a linear wavepacket is local- 
ized in a finite region of space by the Anderson localization 
mechanism. Noise is instead known to break the coherence 
that is necessary to achieve localization, giving rise to a dif- 
fusive expansion of the wavepacket, as predicted by several 
theoretical approaches 1 16-20 1 and also observed in experi- 
ments with atoms and photons GTH251 . Finally, also a weak 
interaction can inhibit the Anderson localization, through the 
coherent coupling of single-particle localized states ll26U36l . 
giving rise to a subdiffusion, i.e. a time-dependent diffusion 
coefficient, that has recently been observed in experiments 
with ultracold atoms [37]. A stronger interaction, in systems 
with limited kinetic energy such as lattices, can instead lead 
to other localization phenomena, such as self trapping l38ll or 
the Mott insulator (39). In this regime the presence of disor- 
der can give rise to new quantum phases B0ll4T1 . However, 
while both microscopic theories and macroscopic models ex- 
ist for noisy or many-body disordered systems, no such theory 



has been developed for the evolution of a wavepacket in dis- 
order with the simultaneous presence of noise and interaction 
effects. Even an intuitive understanding of the problem is pre- 
vented by the difficulty in combining the incoherent dynamics 
generated by noise with the coherent coupling due to interac- 
tion. 

In this work we employ an ultracold Bose-Einstein conden- 
sate in a disordered optical lattice to investigate the general 
features of the expansion of a wavepacket subjected to 
controlled broadband noise and weak repulsive or attractive 
interaction. In general, we find that the combination of noise 
and interaction gives rise to a faster, anomalous diffusion 
of the wavepacket, with a time-exponent that depends in a 
complex manner on the system parameters. Our extensive 
exploration, supported by numerical simulations, indicates 
that such dynamics can be modeled quite accurately over a 
wide range of parameters by a generalized diffusion equation 
in which the two diffusion terms due to noise and interaction 
simply add up. This surprisingly simple result seem to 
persists also in regimes where a perturbative modeling of the 
individual diffusion mechanisms breaks down. For partic- 
ularly strong interaction strengths or narrowband noise, we 
find instead a more complex interplay of noise and interaction 
in the expansion, which cannot be modeled in a simple 
way and will require a more detailed investigation. Besides 
providing a general model for the diffusion in an interacting 
noisy system, this study highlights the capability of quantum 
gases in disordered optical potentials to investigate other 
general problems related to noise in quantum systems. 

Results 

Dynamics in a disordered lattice. The experiment is 
based on a Bose-Einstein condensate of 39 K atoms in the fun- 
damental energy band of a quasiperiodic potential, which is 
generated by perturbing a strong primary optical lattice with 
a weak, incommensurate, secondary one (Fig. [lj). The site- 
to-site tunneling energy J of the main lattice and the disorder 
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FIG. 1. Expansion with disorder, noise and interaction, a) Noisy 
quasiperiodic potential (solid line), realized by a static main optical 
lattice overlapped to a secondary one (dashed line), with a broadband 
amplitude modulation, b) Measured time evolution of the width for 
different values of noise amplitude and interaction energy for A/ J = 
4: A — and E int ~ (stars), A = 0.9 and E int ~ (triangles), 
A = and E int ~ 0.8J (squares) and A = 0.8 and £ ini ~ 0.8J 
(circles). The lines are fits with Eq. {2j. 



strength A characterize the corresponding Hamiltonian 
In the non-interacting case the system shows an Anderson- 
like localization transition for a finite value of the disorder 
A = 2 J. Above this threshold all eigenstates are exponen- 
tially localized, with a localization length £ w e?/ln(A/2J) 
(43), where d is the spacing of the main lattice. 

The noise is introduced by an amplitude modulation of 
the secondary lattice, with controllable strength A and with 
a frequency uj m that is randomly varied in a proper interval 
(Methods). This corresponds to a broadband spectrum with a 
controllable width of the same order of the energy bandwidth 
of the quasiperiodic potential, i.e. W ~ 2A + 4 J. We note 
that our noise is non-dissipative, and in general we have 
an out-of equilibrium situation in which the fluctuation- 
dissipation relation does not hold. However, the finite 
bandwidth of the lattice sets a limit of the order of 4 J to the 
maximum kinetic energy that can be pumped into the system 



by the noise source. Numerical simulations we performed 
indicate that after a typical time t — 0.1 — Is our system 
has reached this limit, and most of the long-time dynamics 
happens in a quasi-equilibrium regime. 

A magnetic Feshbach resonance allows to control the con- 
tact interaction l44l . by tuning the s-wave scattering length a 
and in turn the mean interaction energy per particle Ei nt ph 
2irh 2 an/m, where n is the mean density and m is the parti- 
cles mass. 

To study the dynamics, we initially prepare the condensate 
close to the ground state of the lattice, also in presence of 
a tight axial trap. We then switch off the trap and study the 
expansion dynamics along the lattice, thanks to an additional 
radial confinement. To do so, we detect the axial density pro- 
file n(x, t) with destructive absorption imaging at increasing 
times t, up to t = 10 s. The measurements we present are 
typically referred to as the root mean square width u of the 
sample along the lattice. 

Observation of anomalous diffusion induced by noise 
and interaction. 

A typical time evolution of the width of the system for A > 
2 J in presence of noise and interaction is shown in Fig.[T] The 
disorder-induced localization is broken by interaction or noise 
alone and also by their combination. In all cases, we observe a 
short-time transient that evolves into an asymptotic behavior, 
which is different in the three cases. To analyse the expansion 
we start by recalling that a solution of the diffusion equation 
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is a Gaussian distribution n(x,t) « exp [—x 2 /2a(t) 2 ], with 
da 2 (t)/dt = D. This would imply a time dependence of the 
distribution width as 



a(t) =a (l + t/t ) a , 



(2) 



with the time-exponent a=0.5. For our experimental distribu- 
tions, we find that a solution of a generalized diffusion equa- 
tion of the form of Eq. |2]) with a as a free parameter, can 
be used to fit the rms width of n(x, t) in all cases. Here <tq 
is the initial width, to 1S a f ree parameter that represents the 
crossover time from the short-time dynamics to the asymp- 
totic regime and the exponent a characterizes the expansion 
at long times. In our data we typically see around one decade 
(in time) of asymptotic expansion. In presence of noise alone 
we typically observe normal diffusion, i.e. a = 0.5, as ex- 
pected in case of white noise. The dynamics is instead subd- 
iffusive, i.e. a < 0.5, in presence of a repulsive interaction; 
as we will discuss later, this is essentially due to a reduction 
of the interaction coupling as the system expands. In pres- 
ence of both noise and interaction we typically observe a non 
trivial anomalous diffusion, with expansion exponent a de- 
pending on the relative value of interaction energy Ei nt and 
noise amplitude A. As we discussed above, we are not aware 
of any theory for this combined problem. We therefore model 
the expansion with a generalized diffusion equation for a 2 (t) 
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where the instantaneous diffusion coefficient is the sum of the 
two coefficients of interaction and noise alone, i.e. 



da 2 (t) 
dt 



D mt (t). 



(3) 



Here D rao j se =const is the diffusion coefficient due to noise 
alone, while Di nt (t) is the time-dependent diffusion coeffi- 
cient for interaction alone. A main result of this work is that 
such general diffusion equation is valid in a wide range of 
parameters, as we will show in detail in the following through 
an analysis of the individual diffusion mechanisms and their 
combination. A similar generalized diffusion equation has 
been theoretically predicted for Brownian motion of classical 
interacting particles [45 1. 

Noise-induced normal diffusion. Let us start by exploring 
in detail the effect of noise on the dynamics of the linear sys- 
tem, which we realize by tuning the scattering length a close 
to zero. As shown in Fig. we observe that a finite noise 
amplitude A ^ results in a slow expansion of the initially 
localized sample. The shape of n(x) keeps being Gaussian 
at all times. We fit a(t) with Eq. from these and other 
data we measure a = 0.45(5), which is therefore consistent 
with normal diffusion, and we extract a diffusion coefficient 
D = a 2 /t . 

We have performed an extensive investigation of the diffu- 
sion dynamics, by exploring different values of the noise am- 
plitude A and the disorder strength A/ J, i.e. different local- 
ization lengths. The measured diffusion coefficients D, shown 
in Fig. [2] are in good agreement with numerical simulations in 
terms of a generalized Aubry-Andre model [30, 46], includ- 
ing a dynamical disorder analogous to the experimental one 
(Methods). 

We can interpret the observed diffusion as an incoherent 
hopping between localized states driven by the broadband 
noise. A perturbative approach suggests D oc T£ 2 and 
r oc A 2 , where the localization length £ represents the natural 
length-scale of the hopping and T is the perturbative transi- 
tion rate |[T9l . For the Aubry-Andre model, in the limit of a 
noise bandwidth equal to the lattice bandwidth, we calculate 
a diffusion coefficient: 



D oc 



A 2 J{Z + df 



h l 



(4) 



which is in good agreement with both the experiment and nu- 
merical simulations (Methods). As shown in Fig. [2] we ob- 
serve the linear dependence of D on A 2 and also an increase 
of D with £. 

The described behavior persists in most of the range of val- 
ues of A and £ we have been able to explore in the experiment. 
The simulations however give a clear indication that large 
values of A and/or £ would bring the system into a different 
regime, where the dependence on the localization length 
becomes weaker. This is expected since, in the presence of a 
sufficiently strong noise, the perturbative approach based on 
localized states must fail. Interestingly, the experimental data, 
the numerical ones and the perturbative model indicate that 
the crossover to this second regime happens at strong noise 
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FIG. 2. Noise-induced normal diffusion: £ and A dependence 
of the diffusion coefficients, a) Measured diffusion coefficient D 
vs the square of the noise amplitude A, for different values of the 
disorder strength. From top to bottom: A/J = 2.5 (rhombuses), 
3 (spheres), 3.5 (stars), 4 (circles), 5.5 (triangles) and 8 (squares). 
The open symbols stand for the experimental measurements and the 
filled ones for numerical simulations. The dashed line represents 
the critical values for the noise strength A c and D above which the 
perturbative approach is no longer valid (Methods). 



amplitudes, A 
lengths, i.e. £ p 



1, for our range of rather short localization 



Interaction-induced subdiffusion. We now discuss the ef- 
fect of the interaction, which is introduced in our system by 
changing the scattering length a to a finite value, hence in- 
troducing a finite Ei nt . The effect of such non-linearity on 
the dynamics on a lattice with static disorder has already been 
studied in theory |26l[32"l and experiments [25 37]. Basically, 
the finite interaction energy breaks the orthogonality of the 
single-particle localized states, weakening the localization. In 
this case one can describe perturbatively the resulting dynam- 
ics as an interaction-assisted coherent hopping between local- 
ized states, with a coupling strength that decreases as the sam- 
ple expands, since the density and hence Ei nt decrease [37]. 
This results in a subdiffusive behavior, i.e. in a diffusion with 
a decreasing instantaneous diffusion coefficient 



D int (t) = 2a i 



"in* 

( — — a{t) 
to 



(5) 



where ai nt is the time-exponent. From da(t) 2 /dt = 
Di n t(t) one gets indeed an evolution of the form 
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FIG. 3. Noise and interaction additive anomalous diffusion, a) Measured time evolution of the width for noise alone (triangles), interaction 
alone (squares), or both (circles), for A/ J = 4. The noise amplitude is A = 0.6, while the interaction strength is Ei nt ~ 0.8 J. The 
experimental data are fitted with Eq. |2j (solid lines). The dashed line is the numerical solution of Eq. {3} using the extracted diffusion 
coefficients from the fits, with the confidence interval shown as a grey area, b) Relative residuals of the fit, for the case with noise and 
interaction, and the solution of Eq. l|3j, with the associated confidence interval (grey area), c) Numerical simulation of the time evolution of 
the width for interaction alone (i), noise alone (ii), and both (iii), for A/J = 2.5, J = 180 Hz, A = 0.195, T — 1 ms and Ei„t = 1.6J. 
The numerical data are fitted with Eq. |2]) (solid lines), whereas the dashed line is the numerical solution of Eq. d) Relative residuals 
of the numerical cr(t) with noise and interaction (iii) and the solution of Eq. The grey area represents the confidence level of 2 standard 
deviations. 



a(t) = <7rj(l + t/to) aint . The precise value of such 
exponent depends on the details of the interaction strength 
and of the spatial correlations of the disorder: for an uncorre- 
lated random disorder and E int ps A a perturbative approach 
predicts a% n t — 0.25. In our experiment we find exponents 
in the range 0.2< a jnt <0.35. A detailed description of the 
various regimes achievable in a quasiperiodic lattice can be 
found in Ref. B71 . 

Generalized diffusion model. When we introduce noise 
and interaction at the same time, we observe an expansion 
that is globally faster than for noise or interaction alone, but 
has an exponent 0.3 < a' < 0.5 which is intermediate be- 
tween the two previous ones. This indicates that both dereal- 
ization mechanisms are playing a role in the expansion. One 
example over many of these observations is shown in Fig. [3] 
which reports three characteristic expansion curves for noise 
or interaction alone, or both. In presence of both noise and in- 
teraction, we also find that the distribution does not manifestly 
deviate from a Gaussian during the expansion. To analyse the 
combined dynamics we use the generalized diffusion equation 
of Eq. j3j, where D noise and D int (i) are separately extracted 
by fitting respectively the case with noise alone and interac- 



tion alone with Eq. d2). The numerical solution of the differ- 
ential equation Eq. (Bb is actually in good agreement with the 
experimental data and also with the numerical results of the 
theoretical model, as shown in Fig. [3] 

These examples of the experimental observations and nu- 
merical simulations strongly support the hypothesis of addi- 
tivity of the two derealization mechanisms, at least at first 
order. In the experimental data (see Fig. [3^-b) we typically 
observe a slightly faster diffusion with noise and interaction. 
This difference is due to the axial excitation in presence of 
noise, which in turn excites the radial degrees of freedom in 
presence of interaction [37] and produces the increased expan- 
sion we observe in Fig. [3J5 (see Supplementary Information 
for details). On the contrary, in numerical simulations, which 
neglect the axial-to-radial coupling, there are much smaller 
deviations from the solution of Eq. (p), of the order of a few 
lattice constants (see Fig. [3}i). 

The good agreement between the experimental data and the 
prediction of Eq. (|3]l persists for the whole range of A and 
Ei n t that are accessible in the experiment, with A ranging 
from 0.4 to 1 and Ej nt adjustable up to s» J. Furthermore, 
as we described above (Fig. |2j, in numerical simulations we 
can explore the region of large A and/or £, where the perturba- 
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tive description of the noise effect fails. Interestingly enough, 
although this regime is not perturbative, we find the additivity 
of noise and interaction mechanisms also in this case. 

An investigation over a broad range of parameters in 
theory indicates that the system should behave in a fully 
symmetric way for attractive and repulsive interaction. We 
have experimentally tested this expectation with a sample 
prepared with attractive interaction (a w-100ao with ao being 
the Bohr radius), where the behavior is fully analogous to 
the one in Fig. [3^-b, i.e. the two sources of derealization 
simply add up. Note that systems prepared with a < have 
not a stable state at low energy, but typically occupy the 
whole energy band from the beginning of the expansion. 
This however does not seem to affect the expansion dynamics. 

Therefore, the large region of validity of Eq. Q supports 
the idea that the observed anomalous diffusion is driven by 
the additivity of the two derealization mechanisms, and 
indicates that the incoherent noise-induced hopping between 
localized states does not destroy their coherent coupling due 
to the interaction. Note that the interaction-assisted diffusion 
tends to vanish as the sample expands, so that one should 
expect a long-time crossover to a regime where interaction 
effects are negligible, and the system diffuses normally due 
to noise alone. While we can clearly observe this crossover in 
the simulations, its experimental characterization is prevented 
by the limited observation time. 

Regimes of noise-interaction interplay. 

While all experimental observations are well described by 
the generalized diffusion model discussed above, our numeri- 
cal simulations indicate that deviations might appear when the 
system is subjected to a trapping mechanism. This can happen 
in both cases of noise with a bandwidth 5v much smaller than 
the system one W/h, or of an interaction that is strong enough 
to produce self-trapping [30, 32, 38 1. In both cases, the simu- 
lations show an interplay of noise and interaction which can- 
not be described by Eq. 

In the first regime, i.e. small noise bandwidth {5v <C W/ h), 
the diffusion is much slower than for a broad noise spectrum 
(8v sa W/h), because just a small fraction of the localized 
states can be coupled by the noise. Very interestingly, the ad- 
dition of the interaction in this regime seems to restore a fast 
expansion, as if the additional coupling mechanism between 
localized states provided by the interaction would effectively 
broaden the noise spectrum. One example of this behaviour 
is shown in Fig. |4] where 8v is 10 times smaller than that in 
Fig. [3]>d. One observes a noise-induced diffusion coefficient 
D that is smaller than in the case of broadband noise. Interac- 
tion is able to restore a diffusion comparable with the one ob- 
tained with interaction and broadband noise, therefore much 
faster than the one predicted by the additivity conjecture. 

Note that in this regime the additivity hypothesis fails de- 
spite of the validity of our perturbative approach to describe 
the noise-induced diffusion. This example and the comple- 
mentary case of strong noise which manifests additivity to the 
interaction, although the perturbative model is not valid, in- 
dicate an independence of the additivity conjecture from the 
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FIG. 4. Deviations from additivity with a narrowband noise. Nu- 
merical simulation of the time evolution of the width (a) and relative 
residuals (b) with the same parameters of Fig.[5]>d, but with T = 10 
ms. Deviations from additivity are evident in the discrepancy be- 
tween the curve with noise and interaction (iii) and the solution of 
Eq. ([3} (dashed line). 

validity of the perturbative description of the noise. 

The second regime, i.e. self-trapping, is a general scenario 
for the expansion in lattices in presence of interaction, which 
is reached when Ei nt is too large to be transformed into 
kinetic energy during the expansion 11301 [38), and the bulk 
of the system is effectively localized. An example is shown 
in Fig. [5] where the self-trapping is achieved by introducing 
an attractive interaction on an initial state that is close to the 
single-particle ground state of the system. Here, we observe 
an interplay of noise and interaction resulting in a breakdown 
of the additivity and in a transient superdiffusive behavior, i.e. 
an expansion with a > 0.5. In absence of noise, the central 
part of the distribution, where Ei nt is large, cannot expand, 
and therefore the rms width a increases slower than normal. 
Noise can break this trapping by providing the necessary 
energy to couple the interaction-shifted states to the lattice 
band. This causes the fast expansion at short times, which 
is even faster than normal diffusion. We note that related 
superdiffusion effects have been predicted in linear systems 
with varying energy |48|. An analogous effect might arise 
when the noise is able to release the large kinetic energy 
associated to the self-trapped part of the system, but further 
studies are needed to clarify this mechanism. 

Discussion 

An effective overview of our extensive investigation can be 
achieved by studying the general behavior of the expansion for 
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FIG. 5. Deviations from additivity in the self-trapping regime. 

Numerical simulations of the time evolution of the width for attrac- 
tive interaction alone (i), noise alone (ii), or both (iii), for A/J = 
2.5, J = 180-ffz, A = 0.4, /3 = -35, T = 0.6 ms. An in- 
termediate super-diffusive (a > 0.5) regime seems to be induced 
by self-trapping, while asymptotically one has normal diffusion with 
exponent a ~ 0.5. The (green) dashed-dotted line shows the diffu- 
sive slope, while the superdiffusive one is represented by the (blue) 
dashed line. 



varying A, Ei nt and 5v in numerical simulations. The quan- 
tity that is convenient to study is the so called participation ra- 
tio, PR=(X)i n i)^ 1 ( see Methods for details), which measures 
the number of significantly occupied lattice sites. The specific 
diagram in Fig. |6| which shows the calculated PR after 10s 
of expansion for a particular Gaussian initial distribution, can 
be used to summarize general properties investigated with our 
system. 

For a broadband noise (Fig. [6}t, Sv « W/h), one finds 
again that Anderson localization is suppressed by either noise 
or interaction, with both leading to an expansion of the sys- 
tem. We remark that the underlaying mechanisms for these 
two expansions are quite different. A weak interaction pro- 
duces a coherent coupling of single-particle localized states; 
the resulting transport can be modeled as a coherent hopping 
between such states, with a characteristic subdiffusive behav- 
ior that is determined by the changing density. Noise instead 
drives an incoherent hopping between localized states, which 
result in a diffusive expansion. Once both transport mecha- 
nisms are combined, one sees an enhanced expansion, as high- 
lighted by the increasing width in Fig. [6^ for finite A and Ei nt . 
Our detailed time-dependent analysis indicates that this com- 
bined diffusion is well modeled by the generalized, additive 
diffusion model of Eq. Q. This surprisingly simple results 
seems not to be limited to the region of parameters where the 
noise- or interaction-induced diffusions can be modeled per- 
turbatively. 

The numerical simulations reveal also a region of stronger 
interaction {Ei nt ~ 4 J), where the interaction itself prevents 
expansion, because of the self trapping-mechanism. This is 
visible as a reduction of the PR for increasing Ei nt in Fig. [6^. 
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FIG. 6. Generalized diffusion regimes. The color scale represents 
the number of significantly occupied sites after 10s of our numeri- 
cally simulated expansion, for A = 2. 5 J, J = 180 Hz and a Gaus- 
sian wavepacket initially distributed over « 13 lattice sites, a) De- 
pendence on A and Eint for a broadband noise (8v ss W/h). b) 
Dependence on Eint and 8v, for a fixed noise strength A — 0.4. 



Here, the addition of noise restores a diffusion, which our 
time-dependent analysis has shown to have an anomalous su- 
perdiffusive nature. In this regime, the model of Eq. ([3]) 
clearly breaks down. 

Finally, an interplay of noise and interaction is seen also 
in presence of narrowband noise, as shown in Fig. [6J3. Here 
the expansion is studied for varying inverse noise bandwidths, 
W/h8v, and Ei nt , for a fixed value of the noise strength 
(A = 0.4). In the non interacting regime, the diffusion con- 
stant decreases for increasing W/h8v. The addition of the in- 
teraction recovers a faster expansion, which however cannot 
be described by Eq. ([3]). 

In conclusion, we have experimentally realized an ultracold 
atomic system to investigate the transport of a wave-packet in 
presence of controllable disorder, noise, and interaction. We 
have used it to characterize the noise-induced diffusion and its 
interplay with an interaction-induced subdiffusion. We have 
observed that the complex anomalous diffusion resulting by 
the simultaneous presence of the two processes, can be re- 
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markably explained by a simple generalized diffusion equa- 
tion in a wide range of parameters. Finally, we have numer- 
ically observed also regimes where the presence of trapping 
phenomena produce a more complex interplay between the 
noise- and interaction-induced diffusion mechanisms. It is in- 
teresting to note that a quasiperiodic lattice might allow to 
study also a different regime that has been extensively ana- 
lyzed in theory fl6l[49]|50l , in which in absence of noise all 
the eigenstates are extended, while the noisy potential is disor- 
dered. This regime can be achieved simply by using a weaker 
disorder, below the localization threshold (A < 2 J). Prelim- 
inary numerical simulations we have performed show that in 
this case the noise slows the dynamics from ballistic to diffu- 
sive. 

Our work has showed how the combination of an ultracold 
atomic system with optical lattices can be effectively em- 
ployed to study quantum phenomena related to noise. Future 
experiments in lattices with higher dimensionality or in 
reduced dimensions will allow to explore the regime of strong 
correlations. Here there is a strong interest in understanding 
the behavior of open quantum systems BP . and how quantum 
phase transitions are affected by noise [52 1. An interesting 
alternative to the use of a noisy optical potential is to employ 
a large sample of a second atomic species with a controllable 
thermal distribution, which is coupled via resonant elastic 
scattering to the atomic system under study. This approach 
might allow to introduce a noise in thermal equilibrium, and 
to simulate phonon-related phenomena ll53l l54l . 
Furthermore, our scheme could be easily used to control the 
temperature of our weakly interacting bosons in presence 
of disorder, to investigate the many-body metal-insulator 
transitions at finite temperatures lf55l . 

Methods 

Quasiperiodic potential. The one-dimensional quasi- 
periodic potential is created by a primary optical lattice com- 
bined with a weaker incommensurate one: 

V(x) = Vi sin 2 (k 1 x + <p x ) + V 2 j 2 sir^ij^x + (p 2 ) ■ (6) 

Here fc, = are the wavevectors of the lattices 

(Ai = 1064.4 nm and A2 = 859.6 nm) and 7 = j± measures 
the commensurability of the lattices. The relevant parameters 
are the spacing d=Ai/2, the tunneling energy J (typically 
150 Hz) of the primary lattice, and the disorder strength A, 
which scales linearly with V 2 |46|. Non-interacting particles 
in the fundamental band of this lattice are described by the 
Aubry-Andre model ll43l which shows a metal-insulator 
transition for a finite value of the disorder A = 2 J [46 1 . 

Interaction energy. The scattering length a is changed 
by means of a broad Feshbach resonance to values ranging 
from a ~ O.lao to about a = 300ao 04). We can define a 
mean interaction energy per atom Ei nt = 2irh 2 an/m, where 
n is the mean on-site density. The atomic sample is radially 
trapped with a frequency uj r = 2ir x 50 Hz. The presence 
of the radial degrees of freedom limits E int to values that 
are not much larger than the kinetic energy J or the disorder 



energy A. 

Noise implementation. In order to have a controllable 
amount of noise, we introduce a time dependence on the am- 
plitude of the secondary lattice potential, which produces a 
time-variation of the on-site energies. In the experiment, the 
noise must be broadband enough to couple to as many states as 
possible within the lattice bandwidth (W/H w (2A + 4 J) /%) 
but at the same time one must avoid excitation of the radial 
modes (pj r = 2ir x 50 Hz) as well of the second band of the 
lattice (AEg ap /h «3-5 kHz). This is achieved with a sinu- 
soidal amplitude modulation, i.e. 

V 2 (t) = Vb(l + Asm(oj m t + 4> m )) , (7) 

with the frequency ui m that is randomly varied in a finite in- 
terval oj m E [ujq — Scj, ljo — 5uj] with a time step T, while 
the phase <p m is adjusted in order to preserve the continuity of 
the modulation and the sign of its first derivative, to avoid fre- 
quency components outside the chosen band. The values used 
in the experiment, whose typical spectrum is shown in Fig. [7] 
are: oj Q /2n = 250 Hz, 6u /2n = 50 Hz, T = 5 ms. At long 
times, t 3> T the resulting noise spectrum decays exponen- 
tially outside the band [ojq — 8w, wo + Su], on a frequency 
scale s» 1/T. Note that at short times, i»T, this spectrum 
contains only few components, which might result in an in- 
efficient excitation of the localized states to provide diffusion 
l56l . In the experiment however we do not see such effect. 
We believe that this is due to an additional broadband back- 
ground noise in the relative phase of the two lattices, which 
is due to fluctuations of the position of the retroreflecting mir- 
ror at acoustic frequencies. We measured this phase noise via 
a Michelson interferometer scheme, and we estimate it to be 
equivalent to a broadband amplitude noise on the secondary 
lattice, about 25 dB below the main noise when A = 1. The 
estimated value is consistent with the small diffusion we have 
for A — and with the diffusion we observe for A > and a 
single frequency modulation, i.e. Soj = 0. 

In the numerical simulations we find, however, a reduced 
expansion until t 3> T. For this reason, we have numerically 
tested also a second noise scheme corresponding to a larger 
width of the power spectrum already at short times and is able 
to better reproduce the experimental combination of the exter- 
nal frequency noise and the uncontrollable phase one. In this 
scheme the frequency ui m is kept fixed, while the phase 4> m is 
randomly varied every time step T. Here, the power spectrum 
is essentially a sine function with width Su> = 2ir/T as soon 
as t > T. 

We have checked with extensive simulations that both 
types of noise lead to normal diffusion at long times. In 
the second case, we observe that the diffusion coefficient D 
grows linearly with the bandwidth of the noise as long as such 
bandwidth is smaller than the lattice bandwidth W/ H. Larger 
bandwidths, i.e. very small T, actually lead to a reduction 
of D, since not all the frequencies can effectively produce 
hopping between the localized states. 

The theoretical model. The numerical simulations pre- 
sented in the paper are based on a generalized version of the 
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FIG. 7. Noise spectrum. Typical noise spectral density for T = 5 
ms and A = 1 ((black-solid line). By including the broadband back- 
ground noise (red-dashed line), due to acoustic noise on the lattices' 
mirrors, we estimate the effective experimental spectrum (grey-solid 
line). 



Aubry-Andre model [43 1, including a mean-field interaction 
term 



*Vi(t) = -(ik+iW + + VMM 2 + Mi\ a 1>i ( 8 ) 

where Vj = A/Jsin(27T7j), /3 is the interaction strength and 
ijjj (t) are the coefficients of the wave function in the Wannier 
basis, normalized in such a way that their squared modulus 
corresponds to the atom density on the j-th site of the lattice. 
All energies are measured in units of the next-neighbour tun- 
nelling energy J, while the natural units for time are h/J. 
Note that this model can describe only the lowest energy band 
of the real system. Excited bands are however not populated 
in the experiment. The interaction parameter f3 in the model 
can be connected to the mean interaction energy per particle 
in the real system, which is defined as 



Ej 



2nh 2 fn(r) 2 d 3 r 
a- 



in 



J n(r)d 3 



(9) 



Here n(r) is the mean on-site density distribution, i.e. the 
solution of the interacting Gross-Pitaevskii problem in a 
single well, normalized to an atom number N/n s , where n s is 
the mean number of sites occupied by the atomic distribution. 
The relation between Ei nt and the interaction strength in the 
model is approximately E int « 2Jf3/n s l30ll46l . 

To quantify the localization we consider two quantities: the 
width of the wave packet measured as the square root of the 
second moment of the spatial distribution \ij}j(t)\ 2 , 



1/2 



and the participation ratio PR. 

PR{t) = 



(10) 



(11) 




FIG. 8. Noise-induced normal diffusion: £ dependence of the 
diffusion coefficients. Normalized diffusion coefficient D/A 2 vs. 
the localization length f, for the regime of small noise strength A 
in both experiment (open squares) and numerical simulations (filled 
squares), and for A = 1 in simulations (triangles). The dash-dotted 
line shows the prediction with the perturbative approach. For numer- 
ical data, solid lines are just guides to the eye. 



measuring the number of significantly occupied lattice sites. 
The quantity (J) represents the average over the spatial 
distribution, defined as (j) = ■ j . 



Perturbative modelling the noise-induced diffusion. In 

this section we consider a perturbative approach to describe 
the diffusion induced by noise. The dynamics of the disor- 
dered system in absence of temporal noise is described by 
Eq. dHJ. In the regime of small noise amplitude (A <C 1), 
we can assume that the only effect of noise is to induce hop- 
ping between different localized states of the imperturbed sys- 
tem. The coupling is driven by the frequency component of 
the noise spectrum that is resonant with the energy difference 
between the states. In the regime of £ « d investigated in 
this work, we can restrict our analysis just to the coupling be- 
tween neighbouring states. The hopping rate depends on the 
noise strength and on the overlap between the noise spectrum 
s{u) and the imperturbed system energy distribution P{y): 



A 2 J 



|(z|sin(27r 7 j)|/)| s 



(12) 



where rj = -j^ J s{y)P{y)dv is < 1. To get an analytic pre- 
diction of the effect of the noise, we approximate the energy 
difference between states as a flat spectrum with the same 
width as the lattice band, W/h (2A + AJ)/h. We ap- 
proximate the overlap | (i\ sin(27T7j) |/) | 2 with | (i\f) | 2 , where 
the states \i) and |/) are assumed to be exponentially lo- 
calized over a distance £. It is also reasonable to consider 
A/J ~ 2e d ^. Finally, when the noise spectrum is a sine 
function with T equal to the inverse of the bandwidth W/ h, 
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we obtain an analytical estimation of the diffusion coefficient 

A 2 J(£ + d) 2 



D = 



m 1 



(13) 



This expression is evaluated for the optimal case of 
T = h/W, which determines the specific numerical prefactor 
1 /3. In the general case we find a more complex dependence 
of D on £, and a coefficient D that decreases with the width of 
the noise spectrum 8v — l/T, i.e. the overlap between s{v) 
and P(y) decreases. As shown in Fig. [8] the diffusion coeffi- 
cient obtained with this heuristic model does almost perfectly 
capture the evolution of D with both A and £. Nevertheless, 
we numerically observe that the perturbation approach does 
not hold for large values of A and £, i.e. £ ^ 2d. In these 
regimes, in fact, the £ dependence of D becomes weaker. 
The good agreement between the perturbative model and 
the observed evolution of D confirms the picture that this 
diffusion is induced by the hopping between localized states, 
in analogy with other disordered noisy systems, such as the 
kicked rotor with noise lfT8ll. 



We can estimate that the perturbative approach fails when 
the energy associated to the perturbation rate becomes com- 
parable with the mean separation energy between states in a 
localization volume, i.e. hT s» A^|. This corresponds to a 
critical noise amplitude 



A r 



1 + e d /« d \ V2 



(1 



|)2 3£ 



(14) 



As shown in Fig.|2 A c and the relative D, calculated respec- 
tively using Eq. ( 14 1 and Eq. ( 13 I, correctly capture the order 
of magnitude for the transition from the perturbative regime 
into a new regime in which the £ dependence of D is weaker. 
On top of that, the fact that A c increases when £ decreases, 
i.e when the disorder is stronger, is also reproduced by our 
simulations. 

There is a further regime of strong disorder that can in 
principle be conceived, in which the localization length is 
smaller than the lattice spacing and d becomes therefore the 
relevant length-scale for the diffusion ifTHl . This is so far not 
accessible in the experiment, because the resulting diffusion 
is too slow to be accurately characterized. 

Supplementary Information 

Noise- and interaction-induced derealization. The con- 
jecture summarized by Eq. (T5J, i.e. that noise- and interaction- 
induced derealization mechanisms cooperate, is found to be 
valid over a broad range of parameters. In the experiment we 
have explored a broad range of values of both A and Ei nt , in- 
cluding attractive interactions. For example, in Fig.|9]we show 
a measurement for an interaction strength similar to the one in 
Fig. [3] but with a larger noise amplitude. The behaviour we 
observe is always the same: the additivity conjecture holds, 
except for a slightly faster expansion in presence of both noise 
and interaction. 




FIG. 9. Noise and interaction additive anomalous diffusion. Mea- 
sured time evolution of the width for noise alone (triangles), interac- 
tion alone (squares), or both (circles), for A/ J = 4 and J = 150 Hz. 
The noise amplitude is A — 0.75, while the interaction strength is 
Eint ~ 0.8 J. The experimental data are fitted with Eq. |2]( (solid 
lines). The dashed line is the numerical solution of Eq. {3}, with the 
confidence interval shown as a grey area. 
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FIG. 10. Increasing kinetic energy in presence of noise. Measured 
time evolution of the radial width ay for E int ~ 0.8J. Each data 
set corresponds to a different value of the noise amplitude: interac- 
tion alone, i.e. A — (squares), A = 0.6 (circles) and A = 0.8 
(triangles). The parameters for the static disorder are A/ J = 4 and 
J = 150 Hz. 



This effect is due to an increased kinetic energy in presence 
of noise, which is coupled to the radial degrees of freedom 
in presence of the interaction [37|. In the experiment we can 
clearly detect the noise-introduced energy by probing the in 
situ radial size of the cloud ay (related to the heating). We 
expect that the average kinetic energy scales approximately 
Fig. 



as at 
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shows that, in the presence of both noise and 
interaction, the heating of the sample is enhanced compared 
to the case of interaction alone. For A = just a small 
heating is present, presumably due to the background phase 
noise. This is due to the fact that the time-dependent potential 
continuously injects energy into the system. 
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